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Abstract 

We have observed the decay B + — > D°K + , using 3.3 million 51? pairs col- 
lected with the CLEO II detector at the Cornell Electron Storage Ring. We 
find the ratio of branching fractions R = B(B+ -» D°K+) /B(B + -» Z) 7r+) = 
0.055 ±0.014 ±0.005. 
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Several authors ]I[ have devised methods for measuring the phase 7 ~ a.rg(V* b ) of the 
Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle, using decays of the type B — > 
DK. Comparison between these measurements and results from other B and K decays 
may be used to test the CKM model of CP violation. CP violation could be manifested in 
B — > DK in the interference between a6->c and a b — > u amplitude (Figure [I]), detected 
when the D meson is observed in a final state accessible to both D° and D°. 

The data used in this analysis were produced in e + e~ annihilations at the Cornell Electron 
Storage Ring (CESR), and collected with the CLEO II detector ||. The data consist of 
3.1 fb -1 taken at the T(4S) resonance, containing approximately 3.3 million BB pairs. To 
study the continuum e + e~ — > qq background, we use 1.6 fb -1 of off-resonance data, taken 
60 MeV below the T(4S) peak. 

CLEO II is a general-purpose solenoidal magnet detector. The momenta of charged 
particles are measured in a tracking system, consisting of a 6-layer straw tube chamber, a 
10-layer precision drift chamber, and a 51-layer main drift chamber, all operating inside a 
1.5 T superconducting solenoid. The main drift chamber also provides measurements of the 
specific ionization, dE/dx, which we use for particle identification. Photons are detected in 
a 7800-CsI crystal electromagnetic calorimeter inside the magnet coil. Muons are identified 
using proportional counters placed at various depths in the magnet return iron. 

We reconstruct D° candidates in the decay modes K + n~ , K + 7c~7T°, or K + 7t~7t + tt~ (ref- 
erence to the charge-conjugate state is implied). Pion and kaon candidate tracks are required 
to originate from the interaction point and satisfy criteria designed to reject spurious tracks. 
Muons are rejected by requiring that the tracks stop in the first five interaction lengths of the 
muon chambers. Electrons are rejected using dE/dx and the ratio of the track momentum 
to the associated calorimeter shower energy. The D° daughter tracks are required to have 
dE/dx consistent with their particle hypothesis to within three standard deviations (a). 
Neutral pion candidates are reconstructed from pairs of isolated calorimeter showers with 
invariant mass within 15 MeV (approximately 2.5a) of the nominal ir° mass. The lateral 
shapes of the showers are required to be consistent with those of photons. We require a 
minimum energy of 30 MeV for showers in the barrel part of the calorimeter, and 50 MeV 
for endcap showers. At least one of the two 7r° showers is required to be in the barrel. The 
7r° candidates are kinematically fitted with the invariant mass constrained to be the 7r° mass. 

The invariant mass of the D° candidate, M(D), is required to be within 60 MeV of the 
nominal D° mass. The M(D) resolution, <Jm(d), is 9 MeV in the K + n~ mode, 13 MeV in the 
K + tt~tt° mode, and 7 MeV in the K + 7t~7t + tt~ mode. The loose M(D) requirement leaves a 
broad sideband to assess the background. 

B + candidates are formed by combining a D° candidate with a "hard" kaon candidate 
track. For each B + candidate, we calculate the beam-constrained mass, M bc = \J E^ — p 2 B , 
where pb is the B + candidate momentum and E^ is the beam energy. M bc peaks at the 
nominal B + mass for signal, with a resolution of cta/ 6c = 2.6 MeV, determined mostly by 
the beam energy spread. We accept candidates with M bc > 5.230 GeV. We define the 
energy difference, AE = E D + \jp\ + M\ — E h , where E D is the measured energy of the 
D° candidate, px is the momentum of the hard kaon candidate, and Mk is the nominal 
kaon mass. Signal events peak around AE = 0, with a resolution of 24 MeV in the K + n~ 
mode, 27 MeV in the K + tt~tt° mode, and 20 MeV in the K + h~'k + 'k~ mode. We require 
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-100 < AE < 200 MeV. 

The largest source of background is the Cabibbo allowed decay B + — > D°it + , distributed 
around AE = 48 MeV. Taking into account correlations between AE and M(D), the AE 
separation between signal and B + — > D°tt + is about 2.3cr in all three modes. The only 
additional variable which provides significant K — tt separation is dE/dx of the hard kaon 
candidate. The dE/dx separation between kaons and pions in the relevant momentum range 
of 2.1 — 2.5 GeV is approximately 1.5a. Our dE/dx variable is chosen such that pions are 
distributed approximately as a zero-centered, unit-r.m.s. Gaussian, and kaons are centered 
around —1.4, with a width of about 0.9. 

Other sources of BB background are B — > D*tt + , B + — > D°p + , and events with a 
misreconstructed D° which pass the selection criteria. Such BB events tend to have low 
AE and broad Mj, c distributions. Continuum e + e~ — > qq events also contribute to the 
background. We reject 69% of the continuum and retain 87% of the signal by requiring 
|cos# s | < 0.9, where 9 S is the angle between the sphericity axis of the B + candidate and 
that of the rest of the event. The sphericity axis, s, of a set of momentum vectors, {pi}, is 
the axis for which I Pi x s l 2 * s minimized. 

In addition to the above variables, discrimination between signal and continuum back- 
ground is obtained from cos8b, where 9b is the angle between the B + candidate momentum 
and the beam axis, and by using a Fisher discriminant The Fisher discriminant is a 
linear combination, T = Yli=i a iViy where the coefficients a, are chosen so as to maximize 
the separation between BB and continuum Monte Carlo samples. The eleven variables, 
are | cos 9 t h r | (the cosine of the angle between the B + candidate thrust axis and the beam 
axis), the ratio of the Fox- Wolfram moments if 2/^0 ; and nine variables measuring the 
scalar sum of the momenta of tracks and showers from the rest of the event in nine, 10° 
angular bins centered about the candidate's thrust axis. Signal events peak around T = 0.4, 
while continuum events peak at T = 2, both with approximately unit r.m.s. 

18.8% of the events have more than one B + candidate, reconstructed in any of the three 
modes, which satisfies the selection criteria. In such events we select the best candidate, 
defined to have the smallest x 2 = [{M bc - M B )/o Mbc f + \{M{D) - M D )/a M (D)} 2 , where M B 
and Md are the nominal B and D masses, respectively. We verify that the distribution of 
the number of candidates per event in the Monte Carlo agrees well with the data. 

The efficiency of signal events to pass all the requirements is 0.4412 ±0.0029 for the K + ir~ 
mode, 0.1688 ±0.0016 for the K+tc-it mode, and 0.2186 ±0.0024 for the K+n-n+n- mode. 
The efficiencies are determined using a detailed GEANT-based Monte Carlo simulation 
and the errors quoted are due to Monte Carlo statistics. 

The number of data events that satisfy the selection criteria, N e , is 1221 in the K + n~ 
mode, 5249 in the K + 7t~tt° mode , and 7353 in the K + 7t~7t + tt~ mode. The fraction of signal 
events in the data samples is found mode-by-mode using an unbinned maximum likelihood 
fit. We define the likelihood function 



e=l 
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(1) 



where Vt(e) is the normalized probability density function (PDF) for events of type t, eval- 
uated on event e, and ft is the fraction of such events in the data sample. The seven event 
types in the sum are 1) signal, 2) B + — > D°tt + , 3) B — > D*ir + ± D°p + , 4) a hard kaon 
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or 5) pion in combinatoric BB events with a misreconstructed D°, and 6) a hard kaon or 
7) pion in continuum events. The fit maximizes £ by varying the seven fractions, ft, subject 
to the constraint J2t ft — 1 • 

The PDF's are analytic, six-dimensional functions of the variables AE, dE/dx of the 
hard kaon candidate, M(D), M bc , T , and cos#b- The PDF's are mostly products of six one- 
dimensional functions, except for correlations between AE, M(D), and M bc in the B + — > 
D°K+ and B + -> D°n + PDF's. 

The dE/dx distributions of K ± ,n ± are parameterized using a Gaussian distribution, 
whose parameters depend linearly on the track momentum. The parameterization is de- 
termined by studying pure samples of kaons and pions in data, tagged in the decay chain 
D* + — > D°7T + , D° — > K~7i + . The parameterization in the other variables is obtained from 
the off-resonance data for the continuum PDF's and from Monte Carlo for the BB PDF's. 

The distribution of B + -> D°K + and B + -> D°n + events in AE - M(D) - M bc space 
is parameterized using the sum of two three-dimensional Gaussians, which are rotated to 
account for correlations. For B — > D*n + + D°p + events we use the sum of two Gaussians 
to parameterize the M bc and AE distributions, and a Gaussian plus a bifurcated Gaus- 
sian for the M(D) distribution. These distributions are essentially uncorrelated due to the 
requirement AE > —100 MeV. For BB events with a misreconstructed D° we use a third- 
order polynomial to parameterize the AE distribution, and a first-order polynomial plus 
a Gaussian for the M(D) distribution. The Gaussian is about three times broader than 
the M(D) resolution, and models the peaking which arises due to the selection of the best 
candidate in the event. The M bc distribution is parameterized using the Argus function || 
f(M bc ) oc M bc Jl — (M bc / E b ) 2 exp[— a(l — (M bc / E b ) 2 )}, plus a Gaussian, which reflects mostly 

B — > D^*\ + or B + — > D°p + events in which we misreconstruct a D°. 

We use a first-order polynomial to parameterize the AE distribution of continuum events, 
and a first-order polynomial plus a Gaussian for their M(D) distribution. The Gaussian 
peaking is due both to real _D°'s and to the selection of the best candidate in the event. The 
M bc distribution is parameterized using an Argus function whose sharp edge is smeared by 
adding a bifurcated Gaussian to account for the beam energy spread. We use the function 

1 — £ cos 2 8b to parameterize the cosOb distributions, and bifurcated Gaussians for the J- 
distributions. 

The results of the maximum likelihood fits are summarized in Table |I[ Averaging over the 
three modes, we find R = B(B + -> D°K + )/B(B+ -> D°tt+) = 0.055 ± 0.014 (statistical). 
This is consistent with the value (fn/fn) 2 tan 2 6 C ~ 0.07, expected from factorization, with 

02 <C a± 0. The x 2 of the average is 1.2 for two degrees of freedom, indicating the consistency 
among the results obtained with the three decay modes. To illustrate the significance of the 
signal yield, contour plots of —2 In £ vs. the number of B + — > D°K + and B + — > D°tt + events 
are shown in Figure 0. The curves represent na contours, corresponding to the increase in 
— 21n£ by n 2 over the minimum value. 

The quality of the fit is illustrated in Figure |3|a, showing projections of the data 
onto dE/dx and AE for events in the B + — > D°K + region, defined by JF < 1.6, 
\M bc - 5280 MeV| < 5 MeV, \M(D) - 1864.5 MeV| < 20 MeV, -50 < AE < 10 MeV, 
dE/dx < —0.75. Requiring that events fall within this B + — > D°K + region reduces the sig- 
nal efficiency by about 50%, but strongly suppresses the background. Overlaid on the data 
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are projections of the fit function. The fit function is the sum of the PDF's, each weighted 
by the number of corresponding events found in the fit and multiplied by the efficiency of the 
corresponding event type to be in the B + — > D°K + region. In Figure |3|b we show projection 
plots for events in the B + — > D°n + region, defined by < AE < 100 MeV, \dE/dx\ < 2.5, 
and with the same requirements on J 7 , and M(D) as in the B + — > D°K + region. These 
projections demonstrate that the fit function agrees well with the data in the regions most 
highly populated by signal and the most pernicious background, and provides confidence in 
our modelling of the tails of the B + — > D°7r + distributions. 

Projections onto Mj, c for events in the signal region (Figure |j) illustrate the relative 
contributions and distributions of signal and background events. Only B + — > D°K + and 
B + — ► D°7r + events peak significantly around Mj, c = Mb, despite the selection of the best 
candidate in the event. 

We conduct several tests to verify the consistency of our result. The fit is run on off- 
resonance data and on Monte Carlo samples containing the expected distribution of back- 
ground events with no signal. In both cases the signal yield is consistent with zero. We 
also fit the data without making use of T or dE/dx, and obtain results consistent with 
those of Table |, with increased errors. We find the branching fraction B(B + — > D°ir + ) = 
(4.82 ± 0.19 ± 0.31) x 10~ 3 , in agreement with previous CLEO measurements ||. The ratio 
between the B — ► D*n + + D°p + and B + — ► D°tt + yields obtained from the fit is consistent 
with the measured branching fractions of these decays fT0" |. In addition, our B + — > D°K + 
result is consistent with that of a simpler, though less sensitive method, used to analyze the 
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Many systematic errors cancel in the ratio R. We assess systematic errors due to our 
limited knowledge of the PDF's by varying all the PDF parameters by ±1 standard deviation 
in the basis in which they are uncorrelated, where the magnitude of a standard deviation is 
determined by the statistics in the data or Monte Carlo sample used to evaluate the PDF 
parameters. The systematic error in R due to Monte Carlo statistics is 0.0033. The error due 
to statistics in the data sample used to parameterize the dE/dx distributions is 0.0028, and 
the error due to statistics in the off-resonance data sample is 0.0017. We assign a systematic 
error of 0.0005 due to the uncertainty in the average beam energy, which we estimate to 
be ±0.16 MeV by using the peak of the Mj, c distribution of B + — ► D°n + events. The total 
systematic error is 0.0047. 



TABLE I. Results of the maximum likelihood fits. Nqk and Nn n are the numbers of 
B + — ► D°K + and B + — > D°ir + events found in the fit, respectively. Errors are statistical only. 
The statistical significance of the signal yield is determined from — 21n£ by fixing the number of 
signal events at zero and refitting the data. 

Mode: K + ir~ K + ir~ir° K + tt~tt + tt^ 

N DK 16.5 ±5.9 13.5 ±8.7 21.5 ± 7.8 

N D7T 240 ± 15 379 ± 22 326 ± 20 

Ndk significance 4.2o" 1.8c 3.8<r 

N DK /N Dn 0.069 ±0.026 0.035 ± 0.023 0.066 ± 0.025 
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In summary, we have observed the decay B + — > D°K + and determined the ratio of 
branching fractions 

B(B + -> D°K + ) 

R = B { {m - DH+ > = 0.055 ± 0.014 ± 0.005. (2) 

Combining this result with the CLEO II measurement B(B + D°tt + ) = (4.67 ± 0.22 ± 
0.40) x 10~ 3 , we obtain B(B + -> D°K+) = (0.257 ± 0.065 ± 0.032) x lO" 3 . 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, Research Corporation, the Natural Sciences and 
Engineering Research Council of Canada, the A. P. Sloan Foundation, and the Swiss Na- 
tional Science Foundation. 
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FIG. 1. Dominant Feynman diagrams of B + — > D°K + and B + — > D°K + . 
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FIG. 2. Contour plots of — 21n£ as a function of Ndk and Njy n , the number of B + — > D°K + 
and £? + — > Z)°7r + events found in the fit, respectively. The dashed line marks the 3d contour, (a) 
D° -> i^+vr", (b) D° -> K+tt-tt , (c) Z)° -> ^r+vr^Tr+vr". Note that the iV D7r axis has a suppressed 
zero. 
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FIG. 3. Projections onto the dE/dx and AE axes of the data (points) and fit function (solid 

curves), summed over the three modes. The dashed and dotted curves are the B + — ► D°K + 

and B + — » Z)°7r + contributions to the fit functions, respectively, (a) B + — > D°K + region, (b) 

5 + — > D°tt + region. 
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FIG. 4. Projections onto the Mb c axis of the data (points) and fit function (solid curve) in 
the B + — > D°K + region, summed over the three modes. Also shown are separate background 
contribtions to the fit function: Continuum, continuum plus combinatoric BB events with a mis- 
reconstructed D°, and all backgrounds (including B + — > D°tt + ). 
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